Numerical estimation of the /3-function in 2D systems with spin-orbit coupling 
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We report a numerical study of Anderson localization in a 2D system of non-interacting electrons 
with spin-orbit coupling. We analyze the scaling of the renormalized localization length for the 2D 
SU(2) model and estimate its /3-function over the full range from the localized to the metallic limits. 



I. INTRODUCTION 



II. MODEL AND METHOD 
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As a general rule, all states in a disordered two di- 
mensional (2D) system of non-interacting electrons are 
localized— There are two exceptions. One is the quan- 
tum Hall effect which occurs in 2D systems subject to 
strong perpendicular magnetic fields, where delocalized 
states exist at the center of a Landau levels Another is 
the metallic phase that occurs in 2D systems with sym- 
plectic symmetry, i.e. in systems with time reversal sym- 
metry but in which spin rotation symmetry is broken. 3,4 
The latter of these, which is the subject of this paper, 
is realized when the spin-orbit interaction is significant: 
this interaction breaks spin rotation symmetry but not 
time reversal symmetry. 

While progress has been m&defi& a satisfactory ana- 
lytic theory of the metal-insulator transition in 2D sys- 
tems with symplectic symmetry has yet to be developed. 
In the absence of such a theory, numerical simulation re- 
mains the most useful tool with which to investigate the 
transitioni 4 *^ 

In a recent Physical Review Letter— we estimated the 
critical exponent for this transition. Prior to our work 
varying estimates of the critical exponent had been re- 
ported. The main obstacle to a higher precision estimate 
of the exponent had been corrections to scaling due to ir- 
relevant scaling variables in many of the models analyzed 
numerically. We overcame this difficulty by proposing an 
SU(2) model for which such corrections are much less 
significant. 

In this paper we present a detailed analysis of the scal- 
ing of the quasi-lD localization length in the SU(2) model 
in the metallic, critical and localized regimes. Scaling 
in the critical and localized regimes was not dealt with 
in our Letter. From this we have estimated the renor- 
malization group /3-function for the quasi- ID localization 
length. Wc also present supplementary results for the 
phase diagram and rule out the occurrence of re-entrant 
behavior. 



A. SU(2) model 

We simulated the SU(2) model^ 

H = e ' C t c i<r - R&jloo'CiaCja' W 

where ct(cj CT ) denotes the creation (annihilation) opera- 
tor of an electron at the site i = (x,y) with spin a. The 
random potential is identically and independently dis- 
tributed with uniform probability on [-W/2, W/2]. Hop- 
ping is restricted to nearest neighbors and the hopping 
matrix R(i; j) is distributed randomly and independently 
with uniform probability on the group SU(2) according 
to the group invariant measure. More explicitly 



R(i;j) 



_ e -27;j sin Pi j 



e «7« sm fo. 



(2) 



with Oij and 7^ uniform on [0, 27r), and /3y distributed 
according to 



p(pw = < t md(3 °i p ~ 5 

1 otherwise. 



(3) 



We perform an SU(2) gauge transformation on the model 
for reasons of numerical efficiency (see Appendix |A"|> . 



B. Transfer matrix method 

We have used two different variants of the transfer ma- 
trix method to estimate the localization length A of an 
electron on quasi-lD strips with transverse dimension Lt 
and length L z . We impose periodic boundary conditions 
in the transverse direction. We use quaternion arithmetic 
to perform the transfer matrix calculations. 

The first traditional transfer matrix method 1 ^ esti- 
mates the localization length by simulating a single very 
long sample L z 3> L t . The length of the sample is in- 
creased until a desired precision for A is obtained. 

The second method-^ called here the ensemble transfer 
matrix method simulates an ensemble of samples with a 
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fixed length L z . Here the number of samples is increased 
until a desired precision for the localization length is 
achieved. To ensure that the estimate of the localiza- 
tion length is independent of L Zl a special choice of the 
distribution of the starting vectors in the transfer matrix 
iteration is required. We take a set of ortho-normal vec- 
tors and perform N r transfer matrix iterations on these 
vectors with Gram-Schmidt orthogonalizations. If N r 
is sufficiently large, a stationary distribution of ortho- 
normal vectors is reached. When vectors sampled from 
this stationary distribution are used as starting vectors, 
the estimate of the localization length becomes indepen- 
dent of L z . 

C. Finite size scaling method 



TABLE I: The results of the scaling analysis in the critical 
region. The best fit parameters for energy E — 1, L t = 
8, 16, 32, 64, 96 and W € [5.2, 6.7] are listed. The precision 
of the data A is 0.1%, except for L t = 96 where it is 0.4%. 
Here, and in later tables, Nd is the number of data points, N v 
the number of fitting parameters, and Q the goodness of fit 
probability. 
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Our analysis is based on the assumption that the renor- 
malized localization length A, defined by 



A A 



(4) 



obeys single parameter scaling (SPS) law that is de- 
scribed by the /MunctioniS 



/3(lnA) = 



din A 
dlni, 



(5) 



In the critical regime, where L t <C £, the SPS hypoth- 
esis implies that 



In A = F 



(6) 



Here v is the critical exponent that describes the diver- 
gence of the localization length at the critical point, F is 
a scaling function, and 



4> = iP(E,W) 



(7) 



is a smooth function of disorder and energy that goes to 
zero linearly at the critical point. Equations © and 
are used to fit the results of numerical simulations for 
systems in the critical regime. Once the form of F and 
the critical exponent v are determined, the /3-function is 
calculated by differentiating Eq. ||BJ|. 

When data outside the critical regime are also included 
in the analysis, it is more practical to use a different form 
of the SPS law that expresses A as a function of the ratio 
of the system size to a single relevant length scale £, 



In A = F± 



Li 



(8) 



The subscript distinguishes the scaling function in the 
metallic and localized phases. We follow the convention 
that + indicates the metallic phase and — the localized 
phase. Data for the metallic and localized phases are an- 
alyzed separately. The /3-function can again be obtained 
by differentiating Eq. JSJ) once F± have been determined. 



We accumulated numerical data for the localized, criti- 
cal and metallic regimes, fitted them with the appropriate 
form. The best fit to the data was determined by mini- 
mizing the x 2 statistic and the quality of the fit assessed 
by the goodness of fit probability Q. The precision of the 
results of the fitting procedure were determined using a 
Monte Carlo method^ and expressed as 95% confidence 
intervals. We did not include any corrections to SPS due 
to irrelevant variables^ since these are negligible for the 
SU(2) model. 9 



III. SCALING ANALYSIS OF THE 
RENORMALIZED LOCALIZATION LENGTH 

A. The critical region 

For the critical region we simulated data with a fixed 
energy. When fitting we approximated the function tp by 
a Taylor series truncated at order 



tp = tpiw + ip 2 w z H h W ' 



where 



w 



W c -W 



(9) 



(10) 



Here W c = W c (E) is the critical disorder for the given 
energy. The function F, which for finite L t is a smooth 
function of energy and disorder, was approximated by a 
Taylor series truncated at order no 



F (x) = In A c + x + a 2 x 2 H + . 



(11) 



The coefficients in both Taylor series, the critical disorder 
and the critical exponent are fitting parameters. The 
results of the scaling analysis are given in Table d and 
the best fit is displayed in Fig. ^ 
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FIG. 1: A plot of the best fit to the data for the critical 
regime. All the curves cross at a common point (the critical 
point) . This illustrates that the magnitude of any corrections 
to SPS are smaller than the precision of the data. 

B. The insulating phase 

For very strong disorder the localization length £ is 
short and our data satisfy L t 3> £. In this limit we expect 
that 

A«f (12) 

Lit 

This corresponds to the following limiting value of the 
/3-function 

lim /3(lnA) = -1 (13) 

Supposing that deviations from this limiting value for 
small but finite A are of the form 

/3(mA) = -1 + aA, (14) 

we arrive at the following form for A 

A- x =o+y (A«l). (15) 

We found that this fits data for the strongly localized 
limit well. The fitting parameters are £ = £ (W) at each 
disorder W, and a. The details of the best fit obtained 
with 1)15(1 are tabulated in Table ITT1 

To fit all the data in the localized regime we used a cu- 
bic spline parameterization of the scaling function. We 
set x = ln(L t /£) and y = In A and fit the data with 
y = /_ (x) where /_ is a cubic spline. The values of the 
function /_ at a given set of x values, the derivatives of 



TABLE II: The scaling analysis for the strongly localized 
limit. The ensemble transfer matrix method with N r = 1000 
and L z — 1000 was used. The best fit to data satisfying the 
criterion A < 1/6 obtained in simulations of systems with 



E — 1 and L t 


£ [24, 128] is 
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The precision 
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4.57 ±.02 




N p 


6 


13.0 


4.07 ±.01 




x 2 


10 








Q 


0.9 



1 




-6 -4 -2 2 4 

ln(L,/§ 

FIG. 2: A plot of the best fit to data in the localized phase. 



/_ at the end points of the interval considered, and the 
localization length for each disorder, are fitting param- 
eters. The subroutines "spline" and "splint" in Ref. 
are used to perform the spline interpolation. To fix the 
absolute scale of the localization length in this spline fit, 
we set the localization length at W = 12 to the value 
obtained with Ill5[l in the strongly localized region. The 
results are tabulated in Table ITTll and displayed in Fig O 



C. The metallic phase 

For weak disorder we follow Ref. and argue as fol- 
lows. When A ^> 1 there is negligible localization of the 
electron in the transverse direction on the quasi-lD sys- 
tem. As a result the electron sees an effective random 
potential that is the average of the random potential in 
the transverse direction. In the longitudinal direction the 
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TABLE III: The fit to data for the localized regime: the 
localization length £ at each disorder W and the parameters 
for the cubic spline interpolation. Note that /_ (x) = F- (e x ). 
The data used are for E — 1, L t £ [16, 128]. The precision of 
A ranges from 0.3% to 1.0%. 
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6.3 


4006 ± 200 


M-6.5) 


0.52 ± .04 
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2061 ± 80 


f (-4) 


0.325 ± .004 


6.5 


1187 ±30 


/-(-2) 


-0.015 ± .004 
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528 ± 10 
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/-(2) 
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/-(4) 


-0.96 ± .04 
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7.12 ± .03 
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FIG. 3: A plot of the best fit to data in the strongly metallic 
limit. 



TABLE IV: The best fit to data for the strongly metallic 
limit. The ensemble transfer matrix method with N r — 10000 
and L x = 10000 was used. The best fit to data satisfying the 
criterion A > 4 obtained in simulations of systems with E = 1 
and L t £ [16, 128] is shown. The precision of A ranges from 
0.3% to 1.0%. 
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€ 




0.0 


1.69 ±.07 


b 


4.48 ± .02 


1.0 


1.99 ±.08 
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0.64 ± .01 


2.0 


3.7 ± .1 






2.5 


5.7 ±.2 


N d 


48 


3.0 


10 (fixed) 


N p 


8 


3.5 


19.7 ±.6 


x 2 


44 


4.0 


45 ± 2 


Q 


0.3 



electron is localized with a quasi-lD localization length 
that can be estimated from perturbation theory for a 
strictly ID system. The result is 



A 



1 

W 1 



(16) 



The limiting value of the (3 function that this corresponds 
to is 



lim /3(lnA) = 



(17) 



For large but finite A we speculate that deviations from 
this can be described by an expansion in powers of 1/A. 
Stopping at the first term we have 



/3(lnA) 



c 

A' 



(18) 



This corresponds to a logarithmic increase of A with Lt 



A = 6±cln^ (A > 1) 



(19) 



Data for large A are well fitted by this form. Here b, c and 
the correlation length £ at each disorder W, are fitting 
parameters. Since the absolute scale of £ in the metallic 
phase is arbitrary, we set the correlation length at W = 3 
to £(W = 3) = 10 to fix the scale. This does not affect 
the form of the scaling function or the /3-function. The 
best fit is tabulated in Table llVl and displayed in Fig. [3] 
To fit data for the whole metallic phase we used a 
cubic spline interpolation of the scaling function. We 
set x = In (Zt/£) and y = In A and fitted the data with 
y = /_)_ (x) where /+ is the cubic spline. The correlation 
length £ at each disorder W and the parameters for the 
cubic spline interpolation are the fitting parameters. The 
best fit is tabulated in Table IV! and displayed in Fig. 21 



IV. ESTIMATION OF THE /3-FUNCTION 

For the critical region, we find after differentiating Eq. 
(|6*|) the following expression for the /3-function 

P(F(s)) = -sF'(s) (20) 
v 

For the metallic or localized phases, differentiating © 
we find 



(3(F ± (s)) = sF±(s) 



(21) 



In all cases the /3-function is easily displayed using a para- 
metric plot. 
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FIG. 4: A plot of the best fit to data in the metallic phase. 



TABLE V: The best fit for the metallic phase: the correlation 
length £ at each disorder W and the parameters for the cubic 
spline interpolation. Note that /+ (a;) = F+ (e x ). The data 
used are for E = 1, Lt € [16, 128]. The precision of the data 
A ranges from 0.3% to 1.0%. 
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1.752 ± .003 


3.0 
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2.03 ± .03 
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0.05 ± .02 
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46 ±2 
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0.07 ± .04 
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119 ±5 
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438 ± 30 


N d 


84 


5.3 


1393 ± 100 


N p 


18 


5.5 


4005 ± 400 


x 2 


81 


5.6 


8223 ± 900 


Q 


0.1 



For the strongly localized and strongly metallic limits 
the appropriate forms of the /3-function in terms of the 
fitting paramaters a, b and c have already been given. 

The resulting /3-function is displayed in Fig. [SJ The 
precise form of the /3-function depends not only on the 
universality class but also on the quantity whose scaling 
is analyzed. The /3-function for the renormalized quasi- 
1D localization length discussed here will differ in detail 
from the /3-function for the mean conductance, or the 
mean resistance, or the typical conductance etc4^i^ It 
will also differ in detail from the /3-function found in 
renormalization group analyses of field theories of An- 



FIG. 5: The /3-function for the renormalized localization 
length. The different regions (strongly localized, localized, 
critical, metallic and strongly metallic) are indicated by the 
alternating use of solid and dashed lines. 



derson localizationAiSiiLi& The only common features 
expected to be shared by all /3-functions for particular 
universality class are the existence of a zero, which sig- 
nals the existence of a transition, and the slope at the 
zero, which is related to the critical exponent. 



THE PHASE DIAGRAM OF THE 2D SU(2) 
MODEL 



The preliminary results for the phase diagram pre- 
sented in our previous paper left open the possibility of 
re-entrant behavior similar to that seen for the Anderson 
modeljiS To determine whether or not such behavior oc- 
curs, the data in our previous paper were supplemented 
by simulations with a fixed disorders W — 1 and W = 2 
and varying Fermi energy. The data were fitted as al- 
ready described in Sections III CI and IIII Al the only dif- 
ference being that in Eq. © we set 



w 



E c -E 



(22) 



and determined the critical energy as a function of disor- 
der E c = E c (W). The results are tabulated in Table IVT1 
and displayed in Fig. H3 Re-entrant behavior is clearly 
ruled out. 
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TABLE VI: The details of the simulations and fits used to 
map out the phase digram of the SU(2) model. 
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FIG. 6: Phase diagram of the SU(2) model. The line is a 
cubic spline interpolation. 



estimated the critical exponent v — 2.746±0.009 and the 
/3-function. We also clarified the phase diagram. 

The properties of the metallic phase in the 2D symplec- 
tic universality class are of particular interest. According 
to the single parameter scaling theory, this phase has a 
prefect conductivityi^Sli This is in spite of the system 
being disordered. This conclusion might be avoided if 
there was some breakdown of single parameter scaling in 
the metallic regime. However, we have verified clearly 
in this work that the renormalized localization length A 
does obey the single parameter scaling law in the metallic 
regime. 



APPENDIX A: SU(2) GAUGE 
TRANSFORMATION 



Here we describe the SU(2) gauge transformation men- 
tioned in the text. The Lyapunov exponents are indepen- 
dent of the choice of gauge. 

Taking the x-direction as the longitudinal direction 
and y-direction as the transverse direction, the local 
SU(2) gauge transformation is given by 



C-xyl 



U(x,y) 



C X yl 



(Al) 



where U(x, y) G SU(2) has elements 



U(x,y) = R{x,y\x- l,y)R(x - 1, y;x - 2, y) ■ ■ ■ 

■■■R(2,y;l,y)R(l,y;0,y).(A2) 

After this transformation, the SU(2) model Hamiltonian 
has the same form as Eq. Q but with the hopping matrix 
R replaced with R, where in the x-direction R(x,y;x + 
l,y) is the unit matrix and in the y-direction 



VI. SUMMARY AND DISCUSSION 

In this paper, we analyzed the scaling of the renor- 
malized localization length in the 2D SU(2) model. We 



R(x,y;x,y + 1) = U{x, y)^R{x, y; x, y + l)U(x, y + 1). 

(A3) 

The matrix R(x, y;x,y + 1) is again uniformly and inde- 
pendently distributed on SU(2). 
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